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,  subroutine  BIN  ART  is  an  Implementation  of  tbs  sslf- critical  estimation  procedure 
of  Paulson,  Presser  and  Lawrence  (1983).  The  logistic,  Gaussian  or  Type  I  extras* 
value  distribution  may  be  selected  as  tolerance  distribution.  Estimates  are  expressed 
in  location -scale  fora  on  entry  and  exit,  but  results  nay  be  printed  out  in  regression 
form,  location -scale  fora,  or  in  both  forms  (see  section  3  for  a  description  af  the 
two  paraaetrisstlons) .  >Xt  is  possible  to  hold  location  parameters  constant  during  the 

estimation  procedure.  Estimation  is  accomplished  by  a  Newton -Raphson  method. 

**  V-,  -  .  “  >.  /  t  ■  ’  ■  .  .  nj  '.  X  *V-  rri  j  vi  • 
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lcatton  (FORTRAN! 


SUBROUTINE  BINARY  (N,  IX,  X,  IA,  NPAR,  ISUB,  ISTART,  IDEP,  IDIST,  C, 
RELTOL,  ABSTOL,  MAXTT,  IPRINT,  IFLAG,  BETA,  XLOGL,  ICOV,  COV,  LNEM, 
MEMORY,  IP AULT) 


3.  Description 


The  routine  is  applicable  to  the  following  modeling  situations  For  1  *  i,...,  n 
(n  is  the  sample  size) ,  let  vj  be  the  stress  variable,  aj  a  zero-one  indicator  of 
withstand  or  failure,  and  Xi  *  (x^, . . . ,  >pi)T  a  column  vector  of  covarlates.  A 
constant  is  Incorporated  as  a  covariate  identically  equal  to  unity.  The  case  p*0  is 
possible,  but  should  be  rare. 

The  tolerance  distribution  and  density,  P(vj)  and  f( vj )  respectively,  depend  on 
the  covarlates,  a  scale  parameter  c,  and  a  vector  of  location  parameters  0  * 
(Si,...,  Sp)r  as  follows,  where 


-  lotlttlc.  »(n)  -  jgjjgfo  . 

,<’1’  "  ITT  ’ 

-  Gaussian,  P(?j)  *  S(ui) 

(♦is  ths  standard  Normal  distribution  function), 

*<*i)  *  - - -  «*P(-%uJ)i 

lo\<S~2n 

•  axtraas  valus,  P(»i)  »  1  -  axp( -sxp< uA ) ) , 

f(vi>  *  777  •*p<u1-s*p<ui)). 

In  prlndpls.  It  is  possible  for  a  to  bs  nsgatlvs,  but  this  contradicts  ths  physical 
notion  of  a  strass  uarlabla. 

m  ths  mathod  of  maxlwim  llksllhood,  paramatars  »  ars  astlmatad  by  maximising 
ths  log  llksllhood. 

*«<•>  *  ^tajlogFtvt)  +  (l-at)  logs(v1)), 

whsrs  S  *  i-P.  Ths  sslf -critical  procsdurs  dspsnds  on  a  usar-spaclflad  quantity  c, 
and  sstlmatss  9  by  solving  ths  systsm 

i£i  fC(Vl)[*1  i|  1°*F< vl )  ♦  d-i)  55  logstvi)]  *  o. 

Whan  c  »  o,  maximum  llksllhood  sstlmatss  ars  obtalnsd.  As  c  lncrsasss  from  ssro, 
lncrsaslngly  robust  sstlmators  rssult.  Ths  ssnsltlvlty  of  paramstsr  ssti-natss  to 
dsparturss  from  modal  assumptions  can  bs  sxamlnsd  by  starting  with  c  *  0  and 
rsfittlng  ths  modal  for  lncraaslng  valuas  of  c.  Nsgatlvs  valuas  of  c  ssam  lass  ussful. 


4.  Hu— rlcal  Method 


For  the  purpose  of  esti— tion,  t he  routine  (internally)  repara— terlses  the 
problem  in  a  "regression”  form .  In  thia  parameterization ,  F(  • )  depends  on 

*1  «  avj  +  M'Xi 

instead  of 


The  repar&meterization  offers  two  advantages: 

1)  Computation  of  the  necessary  partial  derivatives  Is  simple; 

2)  For  maximum  llkehood  estimation  with  a  logistic  tolerance  distribution. 

It  can  be  shown  that  the  regression  parameterization  results  in  a  concave 
irtriaHna  problem.  We  anticipate  that  it  will  have  fairly  good  properties 
In  the  —re  complicated  cases. 

It  has  a  disadvantage  In  ter—  of  potential  111 -conditioning  of  the  Hessian  (or 
Jacobian)  matrix,  so  that  the  Input  data  should  be  sensibly  scaled  (see  Section  11) . 
Paras— ers  are  assumed  to  be  expressed  In  the  easier  to  understand  location -scale 
form  on  entry,  and  are  transformed  back  to  location -scale  form  on  exit,  so  the  user 
need  not  worry  about  details  of  reparameterization . 

Let 


Set  «  «<s( ) [at  j|  logF(Vi)  +  (l-at)  logS(vt)]  . 


The  gradient  used  in  Newton -Raphson  iteration  has  0th  component 

n  n  4 

g b  *  n*1  £  Set*  and  the  Hessian  has  (©,  S’)  component  Hee*  ■  n"1  £  — ,  Set, 

i*l  i*i  99 

(the  scaling  by  n*1  should  be  helpful  when  n  is  large). 


When  c  *  o ,  the  estimated  asymptotic  covariance  matrix  of  tha  estimators  la 
n‘1(-H)*1,  tame  whan  c  s  o  It  la  n*lH'lVH*A,  whara  V  has  (9,9*)  componant 
n 

v90*  3  n'1  £  3eise' l •  Tha  aatlmatad  asymptotic  covariance  matrix  la  ax¬ 

ial 

praams  d  in  location -seals  form  by  pra-  and  post  -  multiplying  it  by  tha  Jacobian  of  tha 
undarlylnt  parameter  transformation. 

Newton -Rapbson  lta ration  converges  whan 

max  |  <  ABSTOL 

and 

e<*> 

max  -  <  RELTOL, 

i  o(^-t) 

whara  suparscripta  raprasant  ltaration  numbs  rs  and  ABSTOL  and  RELTOL  ara 
uaar-suppUad.  Tha  uaar  also  spaclflaa  a  maximum  allowmbla  numbar  MAXIT  at 
ltaratlons. 

Solution  of  llnaar  aquations  and  matrix  lnvsrslon  la  accompllahad  by  subroutlnas 
OECOMP  and  SOLVE ,  takan  from  Forsyths,  Malcolm  and  Molar  (1977).  Thasa  routlnas 
ara  of  high  numarical  quality,  and  provida  an  astlmata  of  tha  condition  numbar  of  tha 
input  matrix,  which  Is  oftan  useful.  m  tha  lntarast  of  portabiUty,  itarativa 

lmprovsmant  la  not  usad. 

5.  Paramatars 

5 . 1  input  Paramatars 

N  -  INTEGER 

Sampla  slsa.  Unchanged  on  exit. 

IX  •  INTEGER 

Row  dimension  of  data  matrix  X.  Unchanged  on  exit. 


X  -  real  array  of  DIMENSION  (IX,  q),  whara  q  a  n.  Data  Matrix 

containing  tha  dapandant  (atraaa)  varlabla  and  all  covarlatas. 
Each  mini  oorraaponda  to  ona  obaarvatlon.  If  paramatars  ara 
bald  flxad,  valuaa  of  tha  dapandant  varlabla  will  ba  chan  gad, 
but  thalr  Input  valuaa  raatorad. 

Onchangad  on  axlt. 

IA  -  INTEGER  array  at  DIMENSION  N. 

withstand  or  fallura  la  lndlcatad  for  obaarvatlon  I  according  aa 
IA(I)  -  0  or  IA(I)  a  o. 

Onchangad  on  axlt. 

NPAR  -  INTEGER. 

Nuabar  of  paramatars  in  tha  nodal. 

Onchangad  on  axlt. 

ISOB  -  INTEGER  array  of  DIMENSION  (NPAR). 

Indlcatas  rows  at  data  matrix  X  corraapondlng  to  paramatars 
(tha  acala  paramatar  la  takan  to  correspond  to  tha  dapandant 
varlabla).  For  lnstanca.  If  ISOB(3)  ■  5,  tha  third  paramatar 
oorraaponda  to  tha  fifth  row  at  X. 
onchangad  on  axlt. 

ISTART  -  INTEGER  array  at  DIMENSION  (NPAR). 

Ihdlcatos  tha  status  at  paramatars  In  tha  nodal,  and  whathar  a 
starting  valua  la  to  ba  suppllad. 
if  ISTART(I)  la  aqual  to 

0,  BETA(I)  la  to  ba  astlaatad,  and  Its  input  valua  la  to  ba 
dlaregardad; 

1,  BETA(l)  la  to  ba  astlmatad,  and  its  input  valua  la  to  ba 
uaad  as  starting  valua; 

2,  BETA(l)  la  to  ba  hald  constant  at  its  input  valua. 

(Saa  S action  n . 2  for  rnmmanta  on  starting  valuaa;  tha  acala  at 
paramatar  la  not  allowad  to  ba  hald  constant. ) 

Onchangad  on  axlt. 

IDEP  -  INTEGER. 

Row  of  dapandant  (atraaa)  varlabla  in  X. 

Onchangad  on  axlt. 

IDIST  -  INTEGER. 

Indlcatas  tha  tolaranca  distribution  daslrad. 

If  IDIST  aquala 

1,  logistic  distribution  will  ba  uaad; 

2,  Gaussian  distribution  will  ba  uaad; 

3,  axtrama  valua  distribution  will  ba  uaad. 

Onchangad  on  axlt. 

C  -  DOOBCE  PRECISION 

Usar -suppllad  constant  for  aalf -critical  astlmatlon. 

Onchangad  on  axlt. 


H^,,:rv»r,tV''LT*  '.t  '«%  'j*  -j*  ,’^j,1*^ 


RELTOL  -  DOUBLE  PRECISION. 

Relative  convergence  tolerance  for  Newton -Raphson  Iteration 
(see  Section  4). 
unchanged  on  exit. 

MAXXT  -  INTEGER. 

Maxima  allowable  n unbar  of  Newton -Raphson  Iteration. 
Unchanged  on  exit. 


IPRINT  -  INTEGER. 

Output  unit  n unbar,  if  IPRINT  a  o,  no  output  trill  be  produced. 
If  IPRINT  >  0,  standard  output  trill  be  produced  on  logical 
output  unit  IPRINT. 

Unchanged  on  exit. 


I  FLAG  -  INTEGER. 

Indicator  for  fom  of  output. 

If  IPLAG  <  0,  the  standard  output  suamary  trill  be  in  regression 
fora.  If  IPLAG  a  0,  output  smarles  trill  be  produced  for 
both  regression  and  location -scale  parane ter lxat Ions,  if 
IPLAG  >  o ,  the  standard  output  susmary  trill  be  In  location- 
scale  for*.  (See  Section  11.3). 

Unchanged  on  exit . 


vl 


.V 


3.2  input/output  (and  associated  dlawnsion )  Paraneters. 


BETA  -  DOUBLE  PRECISION  array  Of  DIMENSION  ( NPAR). 

On  entry,  contains  starting  values  as  specified  by  I START. 

On  exit,  contains  parasmter  estljaates,  in  location  scale  for*. 

XLDGL  •  DOUBLE  PRECISION. 

On  exit,  contains  the  log  likelihood  If  c  *  0  (maximal 
likelihood  estimation),  and  zero  otherwise. 

ICOV  -  INTEGER. 

Row  dimen  is  Ion  for  COV. 

Unchanged  on  exit. 

COV  -  DOUBLE  PRECISION  array  of  DIMENSION  (ICOV,  q),  where  q  a  NPAR. 

On  exit,  estimated  asymptotic  covariance  information  for  the 
parameters.  The  dlagnal  contains  standard  errors,  the  strict 
lower  triangle  correlations,  and  the  strict  upper  triangle 
covariances.  If  a  parameter  is  held  constant,  all  correspond¬ 
ing  entries  are  zero,  the  covariance  matrix  will  be  for  the 
location -scale  parameterization  rulers  IPLAG  <  o,  when  It  will 
be  In  regression  form.  See  Section  11.3. 


dimension)  purmtKi 


LMEM  -  INTEGER. 

Length  of  work  array  MEMORY,  aa  declared  In  calling  program  unit. 
LMEN  a  MPAR  +  4*MACT*(  1+MACT)  +  2*max(NFDC,  2*NACT),  Where 
NACT  la  the  nuaber  of  parameters  estimated  and  NFDC  a  NPAR -NACT 
is  the  number  at  parameters  held  constant. 

Unchanged  on  exit. 

MEMORY  -  INTEGER  array  of  DIMENSION  (LMEM) . 

Ueed  as  workspace. 


5.4  Diagnostic  Parameter 
IPAULT  •  INTEGER. 

Unless  the  routine  finds  an  error  or  gives  a  warning,  IPAULT 
will  be  0  on  exit.  See  Section  6. 


6.  Error  indications  and  Warnings 

Errors  or  warnings  specified  by  the  routines 

IPAULT  <  o  If  IPAULT  *  -I,  I  >  o,  then  an  arithmetic  exception  was  about  to 
occur  on  observation  I,  while  computing  partial  derivatives.  This 
failure  Mould  be  rare.  Jt  the  data  have  been  sensibly  scaled,  and 
the  starting  values  are  not  too  bad  (see  sections  11.1-11.2) ,  the 
probable  cause  Is  a  data  error.  The  routine  stops  as  soon  as  the 
error  Is  detected,  and  output  parameter  values  are  not  of  interest. 

IPAULT  »  *  input  parameter  outside  expected  range.  This  failure  will  occur  If, 

on  entry,  N  <  1,  NPAR  <  1,  DC  <  NPAR,  ICOV  <  NPAR,  RELTOL  *  0, 
AB3TOL  <  0,  MAXXT  <  1,  IDI5T  <  1,  IDI3T  >  3,  IPLAG  *  0  and 
IPRINT  a  o  (see  Section  11.3),  ISTART(I)  <  o  or  ISTART(I)  >  2  for 
some  I,  ISTART(I)  ■  2  for  all  I,  ISUB(I)  *  IDEP  for  all  I,  or 
ISUB(I)  ■  IDEP  and  ISTART(I)  a  2  for  some  I.  (The  last  restriction 
means  that  the  scale  parameter  cannot  be  held  constant. )  The 
routine  stops  without  doing  any  calculations. 

IPAULT  >  2  Insufficient  workspace.  This  failure  will  occur  If,  on  entry, 

LMEM  <  NPAR  +  4*NACT*(  1+NACT)  +  2*MAX(  NPDC,2«NACT) ,  where 
NACT  Is  the  number  of  parameters  estimated  (ISTART(I)  <  2),  and 
NPIX  a  NPAR  -  NACT  Is  the  number  of  parameters  held  constant 
(ISTART(I)  «  2).  The  routine  stops  without  doing  any  calculations. 

IPAULT  «  3  The  Besslan  matrix  has  become  numerically  singular  during  Newton  - 
Raphson  iteration.  The  routine  transforms  estimates  to  location - 
scale  form,  restores  values  of  the  dependent  variable  If  parameters 
were  held  fixed,  and  stops,  output  parameter  values  are  not 
generally  of  Interest. 


IPAULT  *  4  The  Newton -Raphson  Iteration  did  not  converge  to  within  RELTOL 
and  ABSTOL  in  the  specified  MAXZT  Iterations .  The  routine 
transferee  estimates  to  location -scale  fora,  restores  values  of 
the  dependent  variable  if  parameter*  were  held  fixed,  and  stops. 
Output  parameter  values  are  not  generally  of  interest. 


7.  Auxiliary  Routines 

SUBROUTINE  PREPAR(  DC,  N,  X,  NPAR ,  BETA,  ISUB,  ISTART,  NACT,  IACT, 
IMAP,  PAR,  NPDC,  IPDCED,  IOEP,  WORK) 

Prepares  for  Iteration  by  setting  up  so**  Indexing  and  working  arrays, 
subtracting  the  effects  of  fixed  parameters  froe  the  dependent  variable, 
and  setting  initial  active  paraaeter  values  in  regression  fore. 

SUBROUTINE  NEWTON(C,  MAXIT,  RELTOL,  ABSTOL,  DIPPER,  XLOGL,  NACT, 
IACT,  PAR,  IPIVOT,  HESS,  HESPAC,  N,  IX,  X,  IA,  IPRINT, 
LWORK,  WORK,  LPAULT) 

Carries  out  the  Newton -Raphson  iteration. 

SUBROUTINE  VARIAB(C,  XLOGL,  REGRES,  DIPPER,  VROUT,  NACT,  IACT, 
IMAP,  PAR,  IPIVOT,  HESS,  HESPAC,  ICOV,  NPAR,  COV, 

DET,  N,  DC,  X,  IA,  WORK) 

Computes  the  estimated  asymptotic  covariance  matrix,  transforms  It  to 
location -scale  form  If  requested,  and  computes  correlations  and  standard 
errors. 

SUBROUTINE  EXPOST( REGRES ,  IDEP,  NPAR,  BETA,  ISTART,  NACT,  PAR,  IMAP, 
N,  DC,  X,  NP EX,  IPDCED,  WORK) 

Sets  the  output  vector  BETA,  transforms  to  location -scale  fora  If  requested. 

If  location -scale  fora  Is  requeted  and  parameters  were  held  fixed,  restore 
Initial  values  of  dependent  variable. 

SUBROUTINE  RESULT  (C,  REGRES,  IPRINT,  IDIST,  N,  IDEP,  NACT,  NPDC, 
XLOGL,  NPAR,  BETA,  ISUB,  ISTART,  ICOV,  COV,  DET) 

Produces  a  standard  output  summary  on  logical  output  unit  IPRINT. 


The  following  routines  are  specific  to  particular  tolerance  distributions,  and  are 
declared  EXTERNAL  in  BINARY.  The  first  three,  prefixed  D,  are  passed  to  NEWTON 
and  VARIAB.  The  next  three,  prefixed  V ,  are  passsd  to  VARIAB. 


SUBROUTINE  DLOGBN(C,  OBJECT,  NACT,  IACT,  PAR,  GRAD,  BESS,  N,  DC,  X, 
IA, ,  ZPAULT) 

SUBROUTINE  DGAUBN  (  -  "  ) 

SUBROUTINE  DEXVBN  (  -  "  ) 

Con  put*  (radiant  and  Haaalan  for  logistic,  Gaussian  and  extras*  value  models, 
reap. 

SUBROUTINE  VLOGBN(C,  NACT,  IACT,  PAR,  V,  N,  EX,  X,  IA) 

SUBROUTINE  VOAUBN( -  -  - ) 

SUBROUTINE  VEXUBN( -  "  - ) 

Compute  v  factor  of  estimated  asymptotic  covariance  matrix  for  logistic, 

Guasslan  and  extreme  value  models,  reap. 

The  following  procedures  perform  general  numerical  tasks.  They  have  been 
Included  In  the  Interest  of  portability,  although  equivalents  exist  an  many  computer 


SUBROUTINE  DECOMP( NDIN ,  N,  A,  UL,  COND,  IPUT,  WORX) 
SUBROUTINE  SOLVE  (NDIN,  N,  A,  B,  X,  IPUT) 


DECONP  decomposes  a  matrix  Into  LU  factors  and  estimates  its  condition,  solve 
solves  a  linear  system,  using  the  results  of  DECONP.  These  routine*  are  In 
Forsythe,  Nalcolm  and  Holer  (1977),  but  the  present  versions  Include  an  extra 
argument  so  that  A  and  B  need  not  be  overwritten. 


SUBROUTINE  DNXNLT(A,  IA,  Nl,  B,  IB,  N2 ,  C,  IC,  N3,  WORX,  LWORX, 
IP  LAG  ,  IPACJLT  ) 


Double  precision  matrix  multiplication  • 

A( Nl  x  N3 )  ■  B(N1  x  N2)  •  C(N2  x  N3),  where  B  is  overwritten  if  IFLAG  <  0 
and  c  la  overwritten  If  IFLAG  >  0. 

DOUBLE  PRECISION  FUNCTION  ALNORN  (X,  UPPER) 

Algorithm  as  66  (Hill,  1973)  to  compute  tall  areas  of  the  standard  Normal 
curve. 

DOUBLE  PRECISION  FUNCTION  RMILLS(X) 

Computes  Z(x)/Q(x),  the  reciprocal  of  Mills'  ratio,  where  x  is  a  standard 
Normal  variate.  This  procedure  is  based  on  procedures  by  Hill  (1973)  and 
Adams  (1969). 


Adam,  A. <3.  (1969) .  Algorithm  39.  A reas  under  the  normal  curve.  Computer  J. . 

12,  197-;196. 

Forsythe,  o.E.  Malcolm,  M.A.  and  Moler,  c.B.  (1977).  computer  Methods  for 
Mathematical  Computations.  Englewood  Cliffs,  Nee  Jersey  1  Prentice 'Ball. 

B1U,  1.0.  (1973).  Algorithm  A3  66.  The  normal  Integral.  Applied  statistics. 

22,  424-427. 

Paulson,  A.S.,  Presser,  M.A.  and  Lawrence,  C.B.  (19*3).  Self -critical,  and 
robust,  procedures  for  the  analysis  of  binary  data.  Operations  Research  and 
Statistics  Report  No.  37-93-Pl,  Department  of  Operations  Research  and  statistics, 
Rensselaer  Polytechnic  institute. 


There  are  no  Internally  declared  arrays. 


10.  Precision.  Machine  Dependent  Constants 


The  routeln  was  developed  on  an  IBM  computers.  The  data  matrix  Is  single 
precision  to  conserve  storage.  TO  convert  to  single  precision,  take  the  following 
steps,  with  the  exception  noted  below  1 


1)  Change  all  DOUBLE  precision  declarations  to  REAL i 

2)  Replace  references  to  double  precision  FORTRAN  library  functions 
with  single  precision  versions,  e.g. ,  EXP  replaces  DEXPt 

3)  Replace  double  precision  constants  by  their  single  precision 
versions,  e.g.,  1.0  replaces  l.ODO. 


Note  1  zn  routines  PREPAR  and  EXPOST,  the  dependent  variable  Is  adjusted  for  the 
effect  of  fixed  parameters.  These  adjustements  must  be  calculated  in  double  precision 


to  avoid  a  loaa  of  significant  digits  in  X  which  could  affect  subssqusnt  computations. 

Procsduros  DLOGBM,  VtOGBW,  DGAUBW,  VGAUBM ,  DEXUBN ,  VEXUBH,  DECOHP, 
ALNORN  and  RMILL3  uss  aachins  -dependent  constants  who—  vlau—  arc  sst  in  DATA 
stats— nts.  Ths—  constants  —  y  havs  to  ba  altsrsd  for  so—  computers.  They  ars 
pointed  out  by  coa— nts  in  ths  progras  units. 

ll.  further  Com— nts 

11.1  Scaling,  conditioning  of  Hessian  — trlx 

The  u— r  should  be  aware  at  potential  proble—  of  ill -conditioning,  narai—  a 
regression -type  para— tarnation  is  used,  calculation  of  the  Hessian  — trlx  involves 
operations  stellar  to  the  formation  of  XTX,  where  X  is  the  data  —trlx. 

Unfortunately,  the  resulting  Hessian  is  often  rather  ill-conditioned.  XT  XPRHTT  >  0, 
an  estl— te  of  the  Hessian's  condition  number  is  printed  out  at  each  Iteration. 
Roughly  speaking,  a  condition  number  in  excess  of  10 T  is  worrieo— ,  although 

condition  numbers  in  excess  of  10iO  can  be  tolerated  when  working  in  double  precision 
on  an  IBM  3081.  The  user  can  and  should  avoid  potentially  excessive  ill -conditioning 
by  scaling  the  data  matrix  X  before  calling  the  routine. 

Although  it  is  not  known  how  to  optimally  scale  a  problem,  it  see—  that  a  nearly 
ideal  scaling  will  be  achieved  if  nonconstant  variables  are  transformed  to  the  range 
C-l,  11,  centered  at  — ro.  However,  such  preci—  scaling  if  often  tedious.  It 

see—  most  important  to  ’’equilibrate"  the  data  matrix  so  that  all  variables  have 
roughly  the  same  (moderate)  magnitude .  such  equilibration  is  often  simply 
accomplished  by  dividing  by  suitable  powers  of  10,  e . g . ,  expressing  voltages  in 

— gavolts  instead  of  kilovolts.  Centering  the  variables  will  further  reduce  tbe 
condition  number.  One  often  centers  covarlates  anyway,  so  that  the  "intercept” 


pmMttr  has  a  claar  Intarpratation . 


11.2.  Starting  valuaa 

3  Inca  Nawton-Raphson  procadura  la  amployad,  starting  valuaa  ara  Important,  if 
no  starting  valuaa  ara  auppllad  by  tha  uaar  (ISTART(I)  *  0  for  all  I),  tha  routlna 
will  uaa  unity  for  tha  acala  paramatar  and  caro  for  all  location  paranatars.  Thaaa 
starting  valuaa  will  not  ba  accaptabla  unlaaa  tha  dapandant  varlabla  haa  baan  cantarad 
and  acalad.  It  la  raconnandad  that  tha  uaar  anploy  tha  naan  and  atandard  davlation 
of  tha  dapandant  varlabla  a a  starting  valuaa  for  tha  “lntarcapt"  paraaatar  (if  any)  and 
acala  paraaatar,  raapactlvaly.  If  tha  nodal  la  aaquantlally  raflt  with  dllforant  valuaa 
of  c.  It  la  racoaaandad  that  aatlmataa  fron  tha  aoat  racant  call  ba  unad  aa  starting 
valuaa  for  tha  naxt  call. 

11.3  Output  flags 

Bacauoa  tha  routlna  will  ganarally  ba  callad  aaquantlally  with  dlffarant  valuaa  of 
c,  it  la  moat  convaniant  to  always  axprass  antry  and  axlt  paramatar  valuaa  in  tha  sama 
forn.  Tha  location -acala  form  la  uaad,  bacauaa  that  form  affords  tha  claarast 
lntarpratatlon .  An  lnconalatancy  arlaas  of  IFLAG  <  o,  for  than  tha  output  astlnataa 
ara  in  location -acala  forn,  but  tha  convarlanca  matrix  la  in  ragraaalon  forn.  Thua, 
If  tha  uaar  wanta  to  naka  aaparata  uaa  of  tha  output  paramatara  of  BINARY,  tha 
routlna  ahould  ba  callad  with  ZPLAG  a  o. 

It  la  rscoamandsd  that  tha  routlna  ba  callad  with  IPRINT  >0.  If  ITLAG  a  o, 
tha  raqulranant  IPRINT  >  0  la  anforcad.  ( whan  iflag  <  0,  tha  raqulranant  la  in 
kaaplng  with  tha  lnconalatancy  mantlonad  abova . )  Tha  standard  output  summary  ahould 
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be  sudden*  for  aoet  applications.  For  the  ■arlaa  flexibility  in  output  and 
Interpretation  of  results,  call  BINARY  with  ZPRXNT  >  0,  iplao  ■  0. 

12.  Exawple 

The  following  prog  ran  illustrates  the  use  of  BINARY,  and  the  standard  output 
produced  when  IFRING  >  o  and  IPLAG  ■  0. 
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'•‘This  paper  describes  the  computational  algorithms  for  computing  generalized 
likelihood  estimators  for  parametric  proportional  hazards  models. 


